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Abstract 
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^ [ We discuss stability of spherically symmetric static solutions in Newto- 

O ' nian limit of Jordan, Brans-Dicke field equations. The behavior of the stable 

Q ■ equilibrium solutions for the spherically symmetric configurations consid- 

er . ered here, it emerges that the more compact a model is, the more stable it 

7' 

u 



is. Moreover, linear stability analysis shows the existence of stable configu- 
rations for any polj^ropic index. 

1 Introduction 



Scalar-tensor theories, in which a long-range scalar field in addition to the 
5^ \ usual tensor fields present in Einstein's theory mediates gravity, have been 

studied in many works (see e.g. [1],[2] ) as natural generalizations of Einstein's 
general relativity (GR). The simplest of them Jordan, Brans-Dicke (JBD) 
theory of gravity [3], [4], in which a scalar field (p acts as the source for the 
gravitational coupling with G ~ was essentially motivated by apparent 
discrepancies between observations and the weak-field predictions of GR. It 
is a well-known fact that most of the mathematical difficulties of theories of 
gravity lie in the high non-linearity of the field equations. However, under the 
special circumstance when the gravitational field is weak one can linearize 
the field equations thereby ignoring this feedback effect. Then the weak-field 
limit is analyzed and the conditions leading to significant deviations of the 
1/r^ Newton's law of gravitation are discussed. 
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The stability of equilibrium configurations for both stars and star clusters 
remains a topic of continuing interest. The fundamental dynamics describ- 
ing galaxies, clusters of galaxies, or globular clusters [5], is a coUisionless 
system. In the general case the coUisionless Boltzmann equation cannot be 
solved because it involves too many independent variables. However, we 
can get certain system with a polytropic state equation, corresponding to 
isotropic velocity dispersion tensors. Moreover, spherically symmetric static 
Newtonian perfect fluid models are the starting point for many discussions 
about stellar structure and evolution. Within this class of models, polytropic 
equations of state have been studied thoroughly [6], [7], [8], [9], [10] 

P = Kp< = Kp^+n' 

where K is non-negative polytropic constant, 7 is the adiabatic index, and 
n the politropic index. In astrophysics pure polytropes have been considered 
and some of the corresponding stellar models have been studied in detail in 
many textbooks [11], although mostly numerically. 

In this paper we consider static spherically symmetric perfect fluid models 
using Newtonian approximations of Jordan, Brans-Dicke (JED) theory of 
gravity. This paper is organized as follows. In section 2 we have constructed 
the field equations for the Newtonian limit of JBD theory. In section 3 
we derive the basic equation with approximation of small linear perturbation 
and give results on stability versus instability for spherically symmetric static 
models. The paper ends with a conclusion in section 4. 

2 Basic equations and their properties 

In this section we discuss the weak-field limit of Jordan, Brans-Dicke theory 
[3], [4]. The JBD theory incorporates the Math principle, which states that 
the phenomenon of inertia must arise from accelerations with respect to the 
general mass distribution of the universe. Differences between predictions of 
JBD theory and observing appeared at study of deciding the field equation: 



2 



Stt 



(2) 



Equations of gravitation field allow significant simplification, if velocity of 
material point far less then velocity of light, so values of order v^/c^ possible 
neglect. In the case of weak field approximation value g^,y must extremely 
little differ from: 

Qixu = 1 for fj,^ ^ 1,2,3, goo ^ -1 
Qnv = for iiy^i^ 

and squares of these deflections possible to neglect. Then 
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Moreover, in static case derived g^i, on time possible to neglect too. Then 
one can change Fgo to oo; or and equations of motion a material point 

(3) for small velocities and weak field takes Newtonian form: 



dU 



where t/ is a gravitation potential, and 



2U 

^00 = 1 



(4) 



Prom equations (1), (2) one can get 
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T, 
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3 + 2a; 



(5) 



For the component 00 equations (5) values of order v/c possible neglect, 
except Too. Component of energy momentum tensor Tqq = pc? consequently 
T = g'^^T^. = g^^Too = -pc\ and 
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1 _ M 2+u 
As far as derived on time and product V^^ we neglect, then 

T3 _ '-'^ 00 
-n-oo — —R — 
ox 

since ^ r^^oo ~ then from (4) 

For scalar potential we have 0;O;o = 2rQQ^ . 
Finely from (1) and (2) 



dz^;(0V[/) = 87r|^p, (6) 



where p is density of mater. Limiting transformation to Newton theory 
of gravitation occurs when | a; | — > oo and (j)— const. 



3 Dynamical stability 

Wc discuss the problem of dynamical stability of the equilibrium solutions 
and consider small time dependent radial perturbations, which still preserve 
spherical symmetry. Let's consider oscillations of configurations consisting 
of ideal gas with adiabatic radial perturbations. Generally the solutions of 
equations of stars oscillations are connected with big mathematical difficul- 
ties, therefore we assume the approximation of small perturbation [12]. Let's 
express w, u and v as radial, latitude and meridianal components of speed in 
spherical coordinates: 
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X — div 




+ 



r sin 9 dip 



1 dv 




then from the equation of continuity 



d_ 
di 



{po + p) 



(Po + p) X, 



(9) 



we have: 




dpo 



dr 



- PoX, 



(10) 



Throughout the paper, we use values with zero designate not perturbed 
values, and without zero a perturbation of the value. In case of ideal gas the 
adiabatic condition will be written down so 



In galactic dynamics 7 is less than 3 [5] which means that no polytropic 
stellar system can be homogeneous. In real stars 7 is the variable and ranges 
from 1 to 00, but possible changes of it are rather small therefore we assume 
for simplicity, that it is a constant. We use for convenience, values e and g 
determined as: 




(11) 




9 = VC/. 



Then from an adiabatic condition we find 




(12) 
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From 



dPo 

and (12) we obtain 

dP 



dr 



Po{gow-e'x). (13) 



dt 

The equations for gravitational potential satisfy 

div {(pogo) = ^^J^^Po 
Perturbation of a gravitational field is given by 



div {(j)go + (pog) = Stt ^^^^^ p, (14) 

The gravitational scalar consist the equilibrium and perturbed part which 
satisfy to the equations 



The equations governing the linear perturbations are obtained by expand- 
ing all functions to first order: 
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Po 



dv 
dt 
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r sin dip 
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[P + Po^) , 



(18) 



dw dP 

Po^^-^-9op-po9- (19) 

Furthermore, we suppose a time dependence of the form 

A{r,t)^A{r) e'"' 

We differentiate (19) on t and (13) on r then taking into account the 
equation of continuity it is received to 



a^w - g'ow - gow' - iag + e^X' + go{l - -f) X = 0, (20) 

where the prime means derivative with respect to r. In a case of radial 
pulsations (8) it transformed to 



X^w'+—, (21) 
r 

The equations (14) and (16) transforms to 



Stt (2 + a;) pow iagocf) 



. 9(f) Snpo , . 

= — T^w. (23) 

9r 3 + 2a; ^ ' 

Then the equation (20) takes a form 



-a^w = -gow H h £ w H 5- + 

00 V r 



+^0 1-7 + — +-w^o + ^^ 24 
V r J r 00 



Performing a change of variable defined by if = r f (r) than from (24) we 
obtain for / 
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e'f'+f (^-7^0)+/ 



'4-37 



9o 



r 00 



3 + 2u; 



(26) 

the pulsation equation transforms to the following Eddington equation in 
the limit a; — > 00 and 6 — const : 



e'f" + f 



- 7^0 + / 



(T^ + (4 - 37) 



90 



0. 



(26) 



Let's put as a first approximation / equal to a constant then from (25) 

'4 — 37 Stt / r pQdr 



90 



+ 



+ 



r (3 + 2u) 4>Q 00/ 



The stability condition one can obtain when put cr^ = 0, in this case 
critical value for 7 is 



_ 4 Stt Jr ppdr 

- 3 + 3(3 + 2a;)0o 30o ' ^ ' 

To begin to understand this problem we have model of gas density dis- 
tribution for clusters of galaxies [5] , with functions of a matter distributions 
for spherically symmetric objects: 



Po = Pc(l+(;^)') ^ (28) 

where the core radius R, central density pc , and the number are parame- 
ters varies from cluster to cluster but has typical value the order of 2/3 which 
implies that the gas mass generally increases linearly with radius. There are 
analogous form of density distribution for globular star cluster [13], in this 
case (3=1. Comparing this relation with a numerical solution of Lane- 
Emden equation for an isothermal sphere [8] one can say that for particular 
case the relative error is less than 5 %. Although these models are simplistic, 
it exhibits many of the key features of the more complex problem. Using 
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relations describing the clusters of galaxies like (28), it is possible to give a 
qualitative answer to the question that reaches critical values of 7. 

Thus, under the assumption of matter distributions (28) we have find 
from (15) expression for not perturbed value of scalar field 0o 



In the previous expressions one can exlain the 0o using the Gauss hyper- 
geometric function. 

For large radius where the effect of the central conditions is very weak 
the solution should asymptotycally approach the exterior solution. It is know 
that the scalar field is a constant outside the matter distributions [14]. In 
empty space there is usual Newtonian universe, but inside "gravitation con- 
stant" depend on matter distributions. Thus, the integration constants Ci 
and C2 arc determined by matching the interior solution (29) to the usual 
exterior Newtonian vacuum solution (j)o=G, 4>q=0. 

To construct a spherically symmetric configurations for given matter dis- 
tributions (28), we choose (3=1 and (3 = 3/2. Knowing 0o and po allows for 
the determination of the critical values of 7 for linear adiabatic radial per- 
turbations of spherically symmetric gas spheres, using equation (27). The 
whole structure of the perturbated gas spheres is thus determined. 

The region of classic Newtonian stability is the rectangle < n < 3. 
However, the conjecture that 7 > 4/3 is a necessary and sufficient condition 
for stability in the Newtonian approximation of Jordan-Brans-Dicke theory 
of gravitation is shown to be false. The critical values of 7 depends on 
the values of the Brans-Dicke parameter uj and on the values of the central 
dencity pc- For both casscs here /3 = 1 and /3 = 3/2 for negative uj is seen 
7 to be less than 4/3 in the instability region. Since scalar field effects are 
stabilizing, one expects stability in wider area than of the Newtonian region. 

4 Conclusion 

In the present work, we have thoroughly analyzed static spherically sym- 
metric configurations in the framework of the Newtonian approximation of 
Jordan-Brans-Dicke theory of gravitation. The classic result for adiabatic 




(29) 



9 



radial pulsations of Newtonian gas spheres is given by the marginal stabil- 
ity 7 > |. Since scalar field effects are stabilizing, one expects stability in 
extended area of the Newtonian region. The behavior of the stable equihb- 
rium solutions for the spherically symmetric configurations considered here, 
it emerges that the more compact a model is, the more stable it is. Further- 
more, in the Jordan, Brans-Dickc theory exist the stable configurations for 
any polytropic index. The similar effect in classic Newtonian theories does 
not exist. Whichever it might be, it is very likely that the same phenom- 
ena could also occur for clusters of galaxies, such as neutron stars and white 
dwarfs. In this sense, the results obtained in this paper can be regarded as 
of a general nature. 

The goal of the present work was to present a stability analyze for model 
of gas density distribution. This configurations can be used to compare with a 
numerical stability code for configurations with polytropic equations of state. 
Together with this, a number of new physical features have been displayed 
concerning the radius-mass relation. 
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